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Effect of weak fluid inertia upon Jeffery orbits 
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We consider the rotation of small neutrally buoyant axisymmetric particles in a viscous steady 
shear flow. When inertial effects are negligible the problem exhibits infinitely many periodic solu¬ 
tions, the ‘Jeffery orbits’. We compute how inertial effects lift their degeneracy by perturbatively 
solving the coupled particle-flow equations. We obtain an equation of motion valid at small shear 
Reynolds numbers, for spheroidal particles with arbitrary aspect ratios. We analyse how the linear 
stability of the ‘log-rolling’ orbit depends on particle shape and find it to be unstable for prolate 
spheroids. This resolves a puzzle in the interpretation of direct numerical simulations of the problem. 
In general both unsteady and non-linear terms in the Navier-Stokes equations are important. 


PACS numbers: 83.lO.Pp,47.15.G-,47.55.Kf,47.10.-g 

Consider a small neutrally buoyant axisymmetric par¬ 
ticle rotating in a steady viscous shear flow. This prob¬ 
lem was solved by Jeffery [T]. He found that the par¬ 
ticle tumbles periodically: it aligns with the flow direc¬ 
tion for a long time and then rapidly changes orientation 
by 180 degrees. There are infinitely many marginally 
stable periodic orbits, the ‘Jeffery orbits’. This degener¬ 
acy means that small perturbations may have substantial 
consequences. It is thus necessary to consider perturba¬ 
tions due to physical effects neglected in Jeffery’s theory. 

For very small particles rotational diffusion must be 
taken into account [2]. The resulting orientational dy¬ 
namics forms the basis for the theoretical understanding 
of the rheology of dilute suspensions [31S]. A second 
important perturbation is breaking of axisymmetry. It 
is known that the rotation of small particles in a sim¬ 
ple shear depends very sensitively on their shape m- 
Third, for larger particles inertial effects must become 
important. This is the question we address here. To 
compute the effect of particle inertia is straightforward 
[SlIS]. But to determine the effect of fluid inertia on the 
tumbling is much more difficult. Despite the significance 
of the question there are few theoretical results, we dis¬ 
cuss them in connection with our results below. 

To understand the effect of fluid inertia on the motion 
of particles suspended in a fluid is a question of funda¬ 
mental importance. But in general it is impractical to 
solve the coupled particle-flow problem, and there is a 
long history of deriving approximate equations of motion 
for the particles, taking into account the unsteady and 
non-linear convective terms in the Navier-Stokes equa¬ 
tions m- The translational motion of a sphere in non- 
uniform flows at low Reynolds numbers, for example, is 
approximately described allowing for unsteadiness of the 
disturbance flow but neglecting convective fluid inertia 
mm- There are many examples where convective fluid 
inertia must be taken into account, leading to drag and 
lift effects mm- In most cases either the unsteady or 
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FIG. 1: Golor online. Spheroid in a simple shear. The flow 
direction is ei, shear direction 62 , and the vorticity points in 
the negative ea-direction. a log-rolling of a prolate particle, 
n is a unit vector along the symmetry axis of the particle, b 
tumbling in the shear plane of a prolate particle, c log-rolling 
of an oblate particle, d tumbling of an oblate particle. 

the non-linear term in the Navier-Stokes equations are 
considered (but see Refs. [T71 [TS]). In our problem both 
unsteady and non-linear convective effects matter. 

We have derived an equation of motion for the orien¬ 
tation of a neutrally buoyant spheroid in a steady shear 
when inertial effects are weak but essential. We show 
how the unsteady and convective terms in the Navier- 
Stokes equations determine the dynamics. Our results 
explain how the degeneracy of the Jeffery orbits is lifted 
by weak inertia. We concentrate on four examples that 
have been discussed in the literature p^H24] : tumbling 
and log rolling of prolate and oblate particles (Fig. [^. 

In this Rapid Communication we give only a brief ac¬ 
count of the formulation of the problem and its pertur¬ 
bative solution (Sections 1 and 2). We focus on the main 
results, Fqs. 0, and 0, and explain their implica¬ 
tions. Details of our calculation are given in Ref. [35]. 



2 


1. Formulation of the problem. Tumbling of a spheroid 
in a simple shear is governed by the shear Reynolds num¬ 
ber Res = safpi/pL (fluid inertia), the Stokes number 
St = (pp//9f)Res (particle inertia), and the particle as¬ 
pect ratio A. Here s denotes the shear rate, pf and Pp are 
fluid- and particle-mass densities, and p is the dynamic 
viscosity of the fluid. We reserve a for the major axis 
length of the particle (used in the definitions of Reg and 
St). The aspect ratio is defined as the ratio of lengths 
along and perpendicular to the symmetry axis. Thus 
A = a/b > 1 (prolate particle), A = 6/a < 1 (oblate 
particle) where 6 is the minor particle-axis length. We 
de-dimensionalise the problem by using the inverse shear 
rate as time scale, particle size a as length scale, and 
ps as pressure scale. For a neutrally buoyant particle 
Res = St. To distinguish the contributions from particle 
and fluid inertia we keep these two parameters separate. 
In dimensionless variables the angular equations of mo¬ 
tion for an axisymmetric particle read 

n — u; A n , (la) 

StL = St(Iuj + iuj)=T. (lb) 

Here n is the unit vector along the particle symmetry 
axis. Dots denote time derivatives, L is the particle an¬ 
gular momentum, I is the moment-of-inertia matrix of 
the particle. The particle angular velocity is lo, and T is 
the torque that the fluid exerts on the particle. To find 
the torque one must solve the Navier-Stokes equations 
for the flow velocitiy u and pressure p subject to no-slip 
boundary conditions on the particle surface 

Res{dtu + u ■ Vu) =—\7p+ y'^u, V-m = 0, (2a) 

u = uj A r ioT r and u = u°° as \r\—>-oo . (2b) 


incompressibility: 

TrS = 0 

symmetry of S: 

=S 

antisymmetry of O: 

= -O 

steady shear: 

00 = -SS, OS = -SO 

normalisation of n: 

n ■ n — 0 

inversion symmetry: 

invariance under n—>■—n, n—>■—fi 


TABLE I: Symmetries constraining the form of Eq. 


taken over the entire fluid vcdume 'f' outside the particle. 
The elements of the matrix U are obtained by solving an 
auxiliary Stokes problem. Details are given in Ref. [23 
Eq. ([^ is exact. The difficulty is that the integrand de¬ 
pends on the sought solution u of Eq. ([^. Therefore we 
follow Refs. [TTI and l20l and evaluate Q to order O(Res), 
the integrand is then only needed to 0(1). More pre¬ 
cisely, we assume that St and Reg are small and of the 
same order, so that ROgSt is negligible. This allows us 
to use the known Reg = St = 0-solutions of ([^ in Eq. ^ . 
The two terms in the integrand in ([^ have the interpre¬ 
tations given in the equation, to linear order in Reg. 

To obtain an equation of motion for n we substitute 
the hydrodynamic torque ([^ into Eq. and expand 

w = -h St -f Reg -h ... . (4) 

Each order in St and Reg must satisfy Eqs. and 
determining the contributions on the r.h.s. of Eq. (4|. 
To lowest order we find the condition = 0. It gives 

= $7 -p An A Sn , (5) 


Here r is a spatial coordinate vector with components 
(ri,r 2 ,r 3 ) in the Cartesian coordinate system 01 , 62,63 
shown in Fig. The undisturbed flow field, u°°, is a 
simple shear flow. We write it as u°° =r 2 ei, so that its 
gradient matrix A has only one non-zero element, Aij = 
Sii6j2. We decompose A into its symmetric part S = 
(A-|-A'^)/2, and its antisymmetric part 0 = (A—A'^)/2. 

2. Perturbation theory. The hydrodynamic torque in 
Eq. ( [Tb| derives from the solutions of Eq. (i). The 
boundary conditions ( 2 b) in turn depend on both par¬ 
ticle orientation n and particle angular velocity u). Thus 
Eqs. ([^ and ([^ are coupled and present a difficult prob¬ 
lem. To proceed we use a reciprocal theorem IIIlEolEe] 
to calculate the torque. Eollowing Ref.ITTlwe find for the 
particular case of a simple shear flow: 


T = -Res dvl]( dtu +u-Vu) 




unsteady 
fluid inertia 


(3) 


The first term in Eq. ^ is the viscous torque com¬ 
puted by Jeffery [T]. The volume integral is the O(Reg)- 
correction to the hydrodynamic torque. The integral is 


where A=(A^ —1)/(A^-|-1) and r2 = (V A u°°)/2, so that 
On = fl An. Eq. ^ is Jeffery’s result [T] for the angular 
velocity of a spheroid in a simple shear, in the absence 
of inertial effects. The second term in Eq. ([^, the St- 
correction, is found to be equivalent to a result given by 
Einarsson et al. j^. We do not reproduce the details 
here because the expression for is lengthy. The 

third term, the 0(Reg)-correction, involves the integral 
in Eq. ([^. But even in perturbation theory [evaluating 
the integrand to order 0(1)] it is difficult to perform the 
integral for arbitrary orientations n. 

3. Symmetries. Exploiting the symmetries of the prob¬ 
lem we can show that it is enough to evaluate the integral 
for only four directions n. The corresponding four inte¬ 
grals suffice to determine the orientational equation of 
motion for n. Here we discuss the idea and give the re¬ 
sulting equation of motion. Details are found in Ref. [25l 

The small-St and -Reg corrections to Jeffery’s equation 
of motion are quadratic in A = 0-|-S. The symmetries 
listed in Table[T|constrain the form of these contributions. 
The resulting equation of motion has only four degrees 
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of freedom which we denote /3i,..., /? 4 : 

n = On+A[Sn—(n • Sn)n] ( 6 ) 

+ /3i(n • Sn)FSn+f32{n ■ Sn)On 
+ /da POSn+/34PSSn. 


The r.h.s. of the first row is Jeffery’s equation, it fol¬ 
lows from Eqs. (la) and ([^. The remaining terms 
are all the terms quadratic in A = O-fS that are al¬ 
lowed by the symmetries listed in Table |l] The pro¬ 
jection P projects out components in the n-direction: 
Pa; = X — (n ■ x)n. The four scalar coefficients / 3 q , are 
linear in St and Reg, and functions of the particle aspect 
ratio: = St/3i®*^^(A)-l-Res/3i^'^'‘^(A). To obtain these 

functions we evaluate Eq. directly for four suitably 
chosen directions n. Comparison with Eq. (§ gives a 
linear system of equations that can be solved for the j3a- 

4 . Results for the coefficients /3a. In two important 
limiting cases the integrand in ([^ simplifies so that we 
can derive explicit formulae for the coefficients Pa- De¬ 
tails are given in Ref. [251 

First, in the limit of large aspect ratios we find that 
particle inertia does not contribute. Pa ' (A) = 0, and we 
obtain that the /3a-coefficients are asymptotic to 


But the correction due to fluid inertia differs from 
the earlier results Eq. (7) in Ref. [191 and Eq. (4.22) in 
Ref. [2T1 In Ref. HH the Navier-Stokes equations ([^ were 
solved iteratively with approximate boundary conditions. 
Only the final result is given, thus we cannot determine 
whether the problem lies in the method or in the algebra. 
We note that Saffman’s assertion that particle inertia can 
be neglected is incorrect, as Eq. ([^ and the results of 
Ref. ini show. We have also verihed Eq. (§ by an inde¬ 
pendent calculation, based on a joint perturbation theory 
in e and Reg using a basis expansion in spherical harmon¬ 
ics. The results are summarised in Ref. [271 and agree with 
Eq. (|^. We also note that Eq. (4.22) of Ref. [2T] violates 
the particle inversion symmetry (Table |T|. 

It follows from Eq. (§ that the unsteady and convec¬ 
tive fluid-inertia terms contribute linearly to This 

enables us to separate their effects to order Reg. For 
large values of the aspect ratio A we find that unsteady 
fluid inertia contributes (81og2A — 12)“^ to / 3 j^®''^ and 
P^<^!>). Comparison with Eq. (j^ shows that the contri¬ 
bution from convective fluid inertia is of the same order. 
For nearly spherical particles, by contrast, we hnd that 
convective inertia dominates (order e), while the contri¬ 
bution from unsteady fluid inertia is smaller, of order . 




for large values of A. The large-A asymptote of Eq. ^ 
agrees with the slender-body limit obtained in Ref. [201 up 
to a factor of Stt. We cannot explain this factor, but have 
verified our results by comparing with an independent 
calculation (Ref. [27], see below). 

Second, we can evaluate the limit of nearly spherical 
particles. We set A = 1/(1 —e) and find to 0(e): 


/3i = 0, /32 = e(St/15 + Reg/35), ( 8 ) 

Ps = e(St/15-37Reg/105), P 4 = e(St/15-bllReg/35). 


In this case particle inertia contributes, and this contri¬ 
bution is consistent with the results of Ref. [21 and also 
with Eqs. (3.15) and (3.16) in Ref. [2T1 


5. Angular dynamics and linear stability analysis. The 
inertial corrections in Eq. ([^ are small in magnitude 
when Reg = St is small, but important because they de¬ 
stroy the degeneracy of the Jeffery orbits. We illustrate 
this effect by analysing four cases: log-rolling along the 
vorticity axis and tumbling in the flow-shear plane, for 
prolate and oblate particles (Fig. [^. In the absence of 
inertial effects these orbits are neutrally stable, as all Jef¬ 
fery orbits in this limit. 

Our analysis is motivated by the fact that recent di¬ 
rect numerical simulation (DNS) results [221 - f21] of the 
problem at small but finite Reg have resulted in a de¬ 
bate as to whether log rolling is stable for prolate parti¬ 
cles, or not. We rewrite Eq. ([^ in spherical coordinates, 
ni = sin0cosy;, n2 = sm9 sinip, 713 = cos0 (the Cartesian 
coordinates are defined in Fig. [^: 


(f= f((p,9) = (A cos 2ip— I)/2 -I- (Pi/ 8 ) sin^ 0 sin 4(/3 — sin 2 ip(P 2 sin^ 9 + P 3 )/A, (9a) 


9 = g{ip, 9) = A sin 9 cos 9 sin (p cos p + sin 


Eqs. (|^ admit two equilibria for 9, log rolling (0 = 0), 
and tumbling in the shear plane (0 = 7 r/ 2 ), see Fig. 

Consider first the linear stability of the tumbling orbit. 
The angle is a monotonously decreasing function of 
time for infinitesimal values of Reg = St. We can thus 


cos 9{Pi sin^ 9 sin^ 2p + P 3 cos 2p + P 4 )/4 . (9b) 


parametrise the orbit by p instead of time, noting that 
p changes from 0 to — 27r during the period time Tp = 
47 r/\/l —A^. We obtain a one-dimensional periodically 
driven dynamical system d9/dp = g{p,9)/f{p,9). We 
define the stability exponent as the rate of separation in 
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one period: 


7t = T ^ lim log |(56>_2^/<56»o| . (10) 

o^o—^0 


Here S0o is a small initial separation from 7r/2 at = 0, 
and 59 -2 -k is the value of this separation at y: = —27r, 
after one period. Linearisation of the 0-dynamics gives: 


" 7o ^d9f{^,7r/2) 


( 11 ) 


Evaluating the integral (111 to order Reg yields an ex¬ 
pression for the exponent 7 t, linear in /3 q,: 


7T = -f + ^-^^(A/l2-^i). (12) 

Log-rolling is a fixed point of the dynamics (§, not a pe¬ 
riodic orbit. But its stability exponent can be calculated 
as outlined above since the (/^-dynamics decouples from 
that of 9, see also Ref. [9]. We find: 


7LR = /34/4. 


(13) 


Using Eqs. 0 and I 
limit (e = (A — 1)/A 


we obtain in the nearly spherical 

0 ) 


7t 

Res 


-2e/21, 


7lr 

Res 


2e/21. 


(14) 


Thus log-rolling is unstable for nearly spherical prolate 
particles (e > 0), and tumbling is stable. For nearly 
spherical oblate particles the stabilities are reversed. An 
earlier approximate theory by Saffman m predicts that 
log-rolling is stable for neutrally buoyant, near-spherical 
prolate spheroids at small Reg, see also Ref.But sta¬ 
ble log rolling has not been observed in DNS for nearly 
spherical prolate spheroids [52H23], and it has been de¬ 
bated how to reconcile this fact with Saffman’s predic¬ 
tion. We have corrected Saffman’s equation of motion. 


As Eq. (14) shows it follows that log-rolling is unstable 


for prolate spheroids at small Reg, consistent with the 
DNS results P5] . 

In the limit of large aspect ratios we find that the ex¬ 
ponents are asymptotic to 


7t 

Reg 


(45 - 30 log 2A)- 


7lr 

Reg 


(15A^)- 


(15) 


We see that tumbling is stable in this limit, log rolling is 
unstable. To determine the stability of the tumbling and 
log-rolling orbits for arbitrary values of A we have com¬ 
puted the /3 q by numerically integrating Eq. ( [^ f or four 
directions n, as outlined above. Figs. 1^ andl^ show 


the resulting exponents. The asymptotes (14) and (15) 


are also shown in Fig. E Fig. Ill demonstrates that log 
rolling is unstable for prolate particles of any aspect ra¬ 
tio. Figs. 0 an also show the separate contributions 
from fluid and particle inertia to the stability exponents. 
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FIG. 2: Color online, a Stability exponent of log-rolling 
(solid line). Separately shown are particle-inertia (dotted) 
and fluid-inertia contributions (dashed). Data computed us¬ 
ing Eq s. (|12|13 1 and numerical results for Pa. (details are given 
in Ref. I25II . Also shown are the asymptotes ( |14| and ( |15[ ), red 
dashed lines, b Same for tumbling in the shear plane. 


We see that the contribution of fluid inertia is in general 
significantly larger than that of particle inertia. 

6 . Concluding remarks. It would be of great interest 
to study by DNS how the stability exponents change as 
Reg is increased, and to determine how the results de¬ 
scribed here connect to those of Ref. [53] valid at larger 
Reg. Second we plan to generalise the calculation sum¬ 
marised here to describe wall effects at small Rcg, by the 
method of reflection [28|. Third, to describe sediment¬ 
ing particles it is necessary to generalise our results to 
Pp 7 ^ pf. Fourth, both the unsteady term and the non¬ 
linear term in the Navier-Stokes equation matter in our 
problem. This raises the question under which circum¬ 
stances both effects matter for the tumbling of small par¬ 
ticles in unsteady flows, and in particular in turbulence. 
Finally we remark that Jeffery orbits are commonly used 
as benchmarks for DNS, despite being valid only in the 
limit Reg = 0. Our solutions provide a new reference 
when fluid inertia is essential but weak. 
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